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Abstract 


The  problem  of  diffraction  of  electromagnetic 
w?. /es  by  an  i-nfinite  slit  in  a  perfect  conductor  at  the 
interface  between  two  different  dielectric  media  is  solved 
by  the  use  of  expansions  of  the  fields  in  each  mediijm  in 
terins  of  appropriate  Mathieu  functions.  3y  iratcliing  these 
expansions  across  the  slit  an  infinite  set  of  linear  equa- 
tions is  obtained  for  the  coefficients  in  the  expansion  of 
the  diffracted  fields.  The  solution  is  put  in  such  a  form 
that  numerical  results  can  be  readily  obtained. 
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!♦  Introduction 


One  of  the  classical  diffraction  problems  is  the  problem  of  diffraction  of 


electromagnetic  waves  by  an  infinitely  long  parallel  slit  in  a  perfectly  conducting 
plane  embedded  in  a  homogeneous  medium.  The  problem  was  first  treated  by  Sieger '-  -^ , 

r2i 

and  extensive  numerical  results  are  due  especially  to  Morse  and  Rubenstein*- -•  and 
Skavlem  L-'-' •  The  general  procedure  is  to  expand  in  terms  of  Mathieu  fimctions  the 
incident,  reflected,  and  diffracted  fields,  and  then  to  evaluate  explicitly  the  coef- 
ficients of  the  series  using  the  boundaiy  conditions. 

In  the  present  paper  we  shall  consider  a  modification  of  this  problem:  we 
assume  that  the  medium  on  one  side  of  the  conducting  plane  differs  in  dielectric  con- 
stant and  magnetic  permeability  from  the  medium  on  the  other  side.  The  general  pro- 
cedure of  expanding  into  a  series  of  Mathieu  functions  is  also  appl-^able  in  this  case, 
but  the  coefficients  in  the  expansion  of  tlie  diffracted  field  can  no  longer  be  evaluated 
explicitly}  they  must  now  be  determined  nijmerically  for  each  wavelength  by  solving  an 
infinite  set  of  linear  equations.  Fortunately  this  infinite  set  reduces,  in  practice, 
to  a  finite  set.  (The  number  of  equations  increases  with  decreasing  wavelength.)  For 
wavelengths  not  smaller  than  one-third  of  the  slit  width,  it  does  not  seem  too  diffi- 
cult to  vjork  out  the  numerical  solution. 
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2,  Formulation  of  the  problem 

Let  z  >  0  be  a  half-space  with  dielectric  constant  e^  and  magnetic  permeabil- 
ity IX.,  and  let  z  <  0  be  a  half -space  with  corresponding  quantities  e,  and  ii^*  At  the 
interface  z  ■  0  there  is  a  perfectly  conducting  plane  with  a  slit»  V^e  introduce  an 
elliptic  cylinder  coordinate  system  by  the  transformation 

(1)   X  =  a  cosh  E,   cos  "ii  ,     y  =  y,       z=a  sinh  4  sin  -9^  ,   (0  <  ^  <oo,  -"  <  ^  <  ")• 

The  slit,  which  is  located  at  -a  <  x  <  a,  is  then  a  singularity  of  the  coordinate 
systems,  because  the  -ralue  of  '♦^  assigned  to  each  point  on  the  slit  may  have  two 
values,* -^J,  r-ienending  on  whether  one  approaches  the  slit  from  the  upper  or  the 
lov/er  half -space. 


Pi.Ture  1 

We  restrict  our  discussion  to  time-harmonic  fields  with  time  dependence  e       , 
where  w  is  the  circular  frequency;  in  general,  we  shall  omit  this  factor  from  our  for- 
mulas.    Maxwell's  equations  then  read 

(2)  icoe  E     =     ctirl  if,  -  ico^  if    =     curl  ^  . 

The  boundary  conditions  are  that 


(3) 

and  that 

(U) 


E     =    E     «=     0         for  z  °  0,      |x|  >  a,     i.e.,  for    '>^  -  0     and    '\  '  +  ^  , 


E  ,  E   ,  H  ,  H     are  continuous  across  the  slit  z  «■  0,    |x|  <  a,  i.e.,  4=0 
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3,  Solution  for  electric  vector  parallel  to  the  slit 

Let  a  plane  wave  be  incident  from  the  direction  designated  in  Figure  1  by 
the  angle  p.  Let  the  electric  vector  of  this  plane  wave  have  only  one  component,  E . , 
which  is  parallel  to  the  slit;  if  we  denote  the  amplitude  of  E.  by  E  then 

(5)  E  =  E  exp  ik,  (x  cos  fi  +  z  sin  j3)|  ,  z  >  0  . 

If  the  slit  were  absent,  we  would  have  a  reflected  wave  with  an  electric 
vector  p>arallel  to  the  incident  wave: 

(6)  E  =  -  E  expMJsLfT  cos  p  -  z  sin  p)  |  ,  z  >  0  . 

Vie  now  expand  each  of  these  waves  into  a  series  of  Mathieu  functions;'^  we 
have  ,  for  0  <  lo  <  n, 

JT  -  2^  i'^ce  (p;  hj)  Mc;l^  (4}  h^)  ce  (^  j  h^) 


(7) 


00 

+ 


and 


(8) 


iii'*l 

F  m 

o        m^O 

^  2r  i^e  (pj  h^)  Ms^^^  (Cj  h^)sej7^}  h^)  , 
ra=l 

where 

(9)  2h^  =  k^a,   k^  «  (^-/e^l"!^   . 

For  the  definitions  and  notations  for  the  ^thieu  functions  used  in  this  paper,  see 
Meixner-Schaef ke ^ -^ ,  For  the  expansion  of  a  plane  wave  see  Section  2,86  of  this  book© 
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To  meet  all  tiie  conditions  of  the  problem  it  will  be  ?een  that  it  is  suffi- 
cient to  assume  that  the  electric  vector  of  the  diffracted  field  also  has  only  a 
single  component,  E,,  which  is  parallel  to  the  slit*  Now  we  expand  also  this  dif- 
fracted field  E,  into  a  series  of  Mathieu  functions: 
a 

(10^    !r  "  ZI  \ti  ^^^    (5J  \^^®m^^  »  ^^         ^  -  °»  ^'^''   °  -  '^  -  "  ' 

o    ni=l 

o    m«=l 

i'he  \xse  of  the  Mathieu  fiinctions  of  the  fourth  kind  arises  from  the  radiation  con- 
dition, nanely  that  the  diffracted  wave  must  be  an  outgoing  wave.  The  omission  of 

terms  with  Mathieu  functions  ce  (''j,  ;  h.  )  and  ce  (''?:  h.,"*  in  (10)  and  (ll)  respective- 

m  J.  m       '     ^ 

ly  is  again  justifie-"  by  the  fact  that  these  terras  are  not  necessary  for  satisfying 
the  reqixLreraents  for  the  solution  of  our  diffraction  problem.     The  parameter  hg  is 
given  by 


(12)  2h2  -     k^a  ,  k^  =    a^/e^T^    . 

The  coefficient  of  the  expansions ,  V.   (i  =1.2)  are  to  be  determined. 

'     xm  '     ' 

On  the  perfectly  conducting  screen  (  '^  °  0,    \  =  ^  ^) »  E.+  E     vardshesj 

therefore  E^  must  also  vanish 'there.     This  is  indeed  the  case  in  e^juations   (10)  and 
d 

(11),  since  the  functions  se     vanish  for    '^    =  0  ai.d   '^   =  i  "• 

In  the  slit  4  =  0  tnt  coriuitiuns   (U)  m\xst  be  satisfied;   that  is,   for  each 
value  of   V    we  must  Pave 


E^    +     E     +     E ,  -    E, 


^  '  o,i<o 


(U) 

i-i-  (E.+  E   +  E,)|  =     l.±E^ 


I  =  0,^<0 


(In  the  above  equations,  the  positive  value  of  '^7    is  used  on  the  left-hand  side  and 
the  sane  value  of  ^    with  a  negative  sign  is  used  on  the  right-hand  side.) 


_  c  - 


We  ncte  that 

^  ■  n  J„  -^    for  S  ■=  C;     se  i-fl)     =  -  se  ('^  )  . 

oz    a  sin  71  04         '      m           in  '• 

Then,  inserting  the  series  expe.rsions  (7),  (8),  (10),  (11)  into  (Ua),  we  get,  for 

00                                 M    -  p                                  CO                                 (.    s                                                                  - 

ZZ,  ^im  K  ^°i  ^^^^m^^  '4^  "    -IZ,  "^2.  ^C'(°'  h^^'J'i  •'  4^  > 
in-±  ni=± 


and 


t^IIiX^Pi   \)  >13^'^'(0;  h^)3e^(^j  h^)  .Z:^  V^«s^^^'(0;  h^)se^(^  ,  h^) 


'J 

These  two  equations  imist  hold  identically  in  "^  •     We  multiply  them  both  by  se  {\\  h- ) 
(n  =  1,2,3,   »••)  and  introduce  the  abbreviations 


Pnm-    l/W^J^^  seJOijh2)d^ 


I'aking  use  of  the  orthonormalization  properties 


/  se^(^i  hj)se^(9^j  h^)d^ 


^  ,  n  -  m 
0  ,  n  /  m 


we  obtain 
(13) 

and 


^in^n^^  (°'  \^     "  -f  ,  Pnm^2m  4^^^°'  ^>  ' 


2  ni=l 


Eliminating  V^  from  (13)  and  (lii),  we  have  finally  the  infinite  set  of  linear  equa- 
tions for  the  V__ 

(15)  mVjPi  h^tef'  CO,  t^)  -C^  p^v^  ^  Ms^^>  (0,.^).  ^tu  („,^)  <'(°"^) 
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From  these  equations  the  V„  can  be  computed  for  given  values  of  h. ,  hp, 
and  S:  and  the  V^  then  follow  from  (13)»  Thus  the  diliracted  wave  is  completely 
determined. 

Finally  we  show  "hat  our  sc3utlon  iTieeti  the  edge  condition,  namely  that 
the  electromagnetic  energy  density  must  be  integrable  over  any  finite  part  of  the 
space.  One  concludes  from  this  that  the  electric  and  the  magnetic  field  strenghts 
may  become  infinite  at  the  edges  z  =  0,  x  =  +  a.  In  the  problem  in  question  we  have 


the  so-called  symmetric  case  of  an  edge,  that  is,  the  angles  around  the  edee,  whjcb 
contain  the  different  media,  are  eq\as.lo  Hence,  as  indicated  by  Karp '-  -'  and  shown  in 

r7i 

detail  by  the  author  in  another  report '--',  the  field  strengths  must  behave  near  the 

edge  as  p   ,  where  p  is  the  distance  from  the  edge  and  t  takes  on  one  of  the  values 

1     3 

«■,  1,  ^, ...  •  Ejrpressing  4  aj^d  i^   near  the  edge  by  p  and  by  the  angle  around  the 

edge,  respectively,  it  can  be  verified  that  in  the  present  problem  the  electric  field 
strength  remains  finite,  since  its  expansion  contains  only  powers  p  ,  p  '   etc., 
while  the  expansion  of  the  magnetic  field  strength  does  in  fact  begin  with  a  term 
containing  p  '  « 
1;,  Magnetic  vector  parallel  to  the  slit 

we  now  consider  the  other  type  of  polarization  that  is,  we  assume  that  the 
magnetic  field  vector  of  the  incident  plane  wave  is  parallel  to  the  slit.  In  this 
case  the  reflected  and  the  diffracted  wave  also  have  their  magnetic  field  vectors 
parallel  to  the  slit.  If  the  amplitude  of  the  incident  wave  is  H  ,  then  the  in- 
cident and  the  reflected  waves  H.  and  H  respectively  are  given  by 

1  r 

H,  ^  -.00 


^    -     axp[ikj_(xcos  p  +  z  sinp)]     -     2^    ^'"ce^^P  J  ^^^V  ^^'  ^^''^ra^  "^  i  ^^ 
o  m«»0 

(16)  ^ 

+  2r    i""  se^(p;  h^)  Ms^^^    (^j  h^)seJX}  i\).  ^  >  0  » 

m^ 


and 
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/    -     exp[ik^(x  cos  p  -  z  sin  p)]     -     2  ^    i™cejp}  h^)Mc^^^e;  hj^)cej^  jhj) 

(17)  ''^         m  2  (l)  9 

-2r"     i'^se^CB;  h,)  Ms^M'^;  i^)se„(9?  I  hf),  z>0 


m=l 


The  difft*acted  wave,  H,,  can  be  expanded  as  follows: 

0  ra«^ 

or 

(19)  j^    «  £3    U^^^""^   (^i  h  )ce  (^  ,  h^),  z  <  0  . 

o  in=0 

2  2 

In  this  e>:pansion  terms  with  fimctions  se  ( i^  j  h- )  and  se  (^  j  h^)  respectively  can 

be  omitted,  as  in  the  preceeding  section.  Consequently,  the  electric  field  component 

parallel  to  the  plane  z  ■  0  vanishes  on  the  screen  C  ■  0,  ^  ■  +  n. 

The  bovmdary  conditions  in  the  slit  are  now 

(hb) 

^  h^\^    "r^  V|c.O,^>0°  '^^"dk  =  0,^<0. 
Inserting  (16),  (17),  (18),  (19)  into  (Ub)  we  obtain 

(20) 

in=U  2  wO 

Now  we  mtiltiply  these  equations  by  ce  (7^  j  h^ )  and  integrate  from  0  to  n  with  respect 
to  '>j^  .  Then  using  the  orthononnaj.zation  relations 
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f     ^^rSl  i  h)^^Jl.  i   ^)d^  = 


for  n  »  m 


0  for  n  /  m 


and  introducing  the  abbreviations 

we  obtain  the  follovaxig  equations: 

(22)    la"ce„(p,  h^)Mo(^)(Oi  h^)  *  Dj^M<:^'''(0,  h^)  -  f    ^^^JK"'   ^°'   "2'  • 

m«0 


.(U) 


.(U)' 


(23)    U^  Mc^-  (Oi  h^)  -  -  -^  n^,  g^U^^  Ifc--  (Oj  h^). 


ew   00 

^2  ri=0 


.(I4)' 


n  -  0,  1,  2,  ., 


By  eliminating  U^  ,  we  arrive  at  the  infinite  set  of  linear  equations  for  the  Up^: 
(2h)  Ui"ce^(pj  Y^)Ko^^hOil^)   ^ Zl^  gJJ^ 


n  ^       En  Mc^"^(0}h,)    /,  N, 
Mc^^')(0-h,)  .-1   ?).  *^   Mci^^  (Ojh^) 


n  ■  0,1,2,  •. .  • 

The  edge  condition  of  the  symnetrical  case  is  satisfied  here  too,  as  is  seen  by  the 

same  argument  as  in  the  previous  case.  The  only  'difference  is  that  now  the  magnetic 

-1/2 
field  strength  is  finite  while  the  electrical  field  strength  behaves  as  p  '   near 

the  edge. 


5«  Nimerical  evaluation 

We  shall  give  here  some  remarks  concerning  the  structure  of  the  enuations 
(l5)  and  the  procedure  for  their  numerical  solution.  These  remarks  apply  ui  essence 
also  to  equation  (23)» 

When  &.  =  Ep,  M-,  "Mp*  ^'^^^   equations  are  greatly  simplified.  In  this  case 
p   "  1  for  n  «  ra,  and  P  =  0  for  n  /  m.  The  equations  (l5)  can  then  be  immediately 


-  9  - 


solved  and  the  diffracted  -wave  coincides  with  the  known  resvilt  for  the  diffraction  by 
a  slit  in  homogeneous  space.  If  cl  /  e^  but  e^n-  •  e-Hp*  "that  is,  if  h-=  h  ,  the 
same  simplification  arises  except  that  the  V^,  are  equal  to  the  V^^  for  the  homogen- 
ecus  case  multiplied  by 


X 


'i*i'2 


C   +  6« 

1    2 


In  the  general  case  e-  (j^  /  ^9^0*   ^^  ^^   convenient  to  introduce  instead  of 


the  V^,  other  unknowns  v  defined  ty 


2m 


m 


m 


Ms^^^'  (Oj  h^)  .  V^^  . 


Then  (l5)  reads 
00 


(25) 


4rrs  nni  m 
m=l 


^2   Ms_^^^  (Oj  hj_) 


n 


Msi>^  (Oj  h^) 


Ms 


ra 
'n 


(0;h^) 


UiV^(PjhJ)Ms^^^'(Oi  h^^)  , 


n  =  1,2,3,  ,,,  , 


One  concludes  from  the  definition  of  the  p   that  if  one  of  the  indices  n.m  is  even 

nin  ' 

and  the  other  is  odd  then  p   ■  0»  So  (25)  separates  into  two  independent  sets,  the 

one  with  n  =  1,3,5,  •••  containing  the  unknoims  v^  ,  v^,  v^^,  ..»,  the  other  with 

n  =  2,U,6,  ...,  containing  the  unknowns  Vp,  v.  ,  V/-,  ...  • 

A  very  rough  estimate  shows  that  the  expression  in  the  brackets  in  (25)  is 

about  of  order  unity,  while  the  order  of  p   is  less  then 

'  nm 


(V2) 


2t 


n-m 


n^tl 


,  with  t  -  -^--^ 


The  right  side  of  (15)  is  of  order  1     »  1 1  \-r')   •  Therefore  the  number  of  equations 
to  be  taken  into  account  in  the  set  (l5)  is  probably  of  order  2h- •  •'•his  estimate  pre- 
supposes that  h^Ap  "   0(1).  For  hn  <  5,  h-  <  5,  we  suppose  that  one  gets  reasonably 
precise  results  with  n  <  10,  that  is,  with  two  sets  of  five  linear  equations,  or  even 
fe>rer« 
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The  numerical  evaluation  of  the  coefficients  in  (25)  can  be  done  conveniently 

M  2 

using  the  'Tables  Relating  to  Mathieu  Functions''-  -'•  There  the  parameter  is  s  =  Uh 

and  certain  quantities  f   ,  g    are  given  as  functions  of  s  from  wMch  we  can  compute 

Ms^^^'(oj  h)  2     r  1 

-^^T-n =  -  i  (g   )^  •  1  -  if 

n    ' 


M^l'Voih)  ./T.    J^   . 


The  p    are  evaluated  by  using  the  expansion 


?     <x>   2t+p   ?    r     1 
^^2t^j>^1i   ^  )  =  11^  Vp  ^^^  3in[(2r+p)^J  ,  P  -  0,1  . 

One  gets 

t2B   2t+p   2^  2n+p  ,   2x 

2t+p 
The  coefficients  B     can  be  derived  from  the  coefficients  Bo^    in  the  above - 
2r+p  '■r-ip 

mentioned  tables;  LJ  we  obtain 


,-^  ,    .„,,.„  [t  ,..„j^-^''^ 


«2^  =   ^°2r.p!A,,  '"-Wj  P'0.1 
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